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Abstract 

In j2] a nonperturbative proof of the g-theorem of Affleck and Ludwig was put 
forward. In this paper we illustrate how the proof of j3] works on the example 
of the 2D Ising model at criticality perturbed by a boundary magnetic field. For 
this model we present explicit computations of all the quantities entering the proof 
including various contact terms. A free massless boson with a boundary mass term 
is considered as a warm-up example. 



*email address: anatolyk@physics.rutgers.edu 



1 Introduction 



A near critical ID quantum system with boundary can be described by a 2D quantum field 
theory defined on an infinite half-cylinder of circumference (3 = 2n/T where T is the system's 
temperature. Let (x, r) : < x < oo, r ~ r + (3 be coordinates on the half-cylinder. We 
put the boundary at x — 0. 



x=0 x—L — oo 



We are interested in systems critical in the bulk but not critical on the boundary. Physi- 
cally the thermodynamic limit is taken by considering a system of finite length L with some 
boundary condition specified at x = L and then sending L to infinity. In the limit L — > oo 
the two boundaries decouple and the leading term in the partition function is 



~ e™ L Wzz' 



where the constants z and z' are the boundary partition functions corresponding to the 
boundary conditions at x = and x = L respectively. If a boundary condition is represented 
by a (properly normalized) boundary state \B) in the quantization in which the Euclidean 
time is along the x direction then z = (B\0) where |0) is the SL(2, (C)-invariant conformal 
vacuum. The partition function Z can be also written as Tr e _/3H ° L where Hql is a Hamil- 
tonian with respect to the Euclidean time chosen along the r direction and the trace is 
taken in the corresponding Hilbert space. Assuming that this Hamiltonian is hermitian Z 
is a nonnegative quantity. Thus for all pairs of boundary conditions the products zz' are 
nonnegartive. Since we are free to multiply all boundary states by a common phase factor 
all z's can be chosen to be nonnegative l . 

We further note that the universal infinite factor e ncL ^ 6/3 and the factor z' drop out from 
normalized correlators of operators inserted in the bulk and/or on the x = boundary. 
Equivalently one can assume that the boundary conditions on a half-cylinder at infinity 
correspond to the bulk CFT conformal vacuum. The bulk stress energy tensor then decreases 
at infinity as 

T^ lk (x,r)^e-^, x^oo. (1) 
The logarithm of the partition function with this boundary condition has the form 

ckL 

In Z = — — + In z L 

6jJ 

and lim^oo zl = z. The boundary entropy is defined as 

s = (l-(3^)\nz. (2) 

Since temperature is the only dimensionful parameter s = s(fij3) where \i is the renormaliza- 
tion energy scale. At a fixed point s = hiz = Ing where g is the universal noninteger gound 

1 This extra care in the discussion of positivity of the boundary partition function z is needed because there is no Tre~@ H 
representation of z. 
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state degeneracy of Affleck and Ludwig lj. It was conjectured in pQ and shown in conformal 
perturbation theory in |2] that s decreases from fixed point to fixed point under RG flow. In 
[3| a nonperturbative proof of this statement which is also known as "g-theorem" was given. 
The proof proceeds via proving a stronger statement - a gradient formula 

-g ab [3 b (3) 



d\ a 



where X a form a complete set of boundary coupling constants, f3 a (\) are the corresponding 
beta functions and the metric g a b is 

P P 

9ab{X) = J Vdn J //rfr 2 (0a(ri)06(r 2 )} c (l - cos[2tt(Yi - r 2 )/0\) . (4) 



Here 4> a ij) are boundary operators conjugated to couplings A a that is 

/idr((j) a (T)) . 



d In z 



o 



Note that throughout the paper we stick to the conventions of |2| in which all operators and 
coupling constants are dimensionless. Hence a factor of [i appears with each integration. 
The gradient formula ©, Q was to a large extent inspired by works on boundary string 
field theory 0j, 0, P, 0. 

In this paper we consider two exactly solvable models: a free boson with a boundary mass 
term and the Ising model at criticality with a boundary magnetic field. The first model was 
first studied in [5] and then used in the study of tachyon condensation in string theory 0. 
It describes a flow from Neumann to Dirichlet boundary condition. The boundary entropy 
for this model was computed in j5] and its monotonic decrease was demonstrated. Because 
of the zero mode the boundary entropy is infinite at the UV fixed point. From this point of 
view this is not a completely clean example of a boundary flow. We consider it in section 3 
as a warm up example, demonstrate the gradient formula and show that the way it works is 
in accordance with the general proof. Our main example is the Ising model with a boundary 
magnetic field. It is considered in section 4. A Lagrangian description of this model was 
given in (111 . It was further studied in ^2] > [S] • in [Ej a local magnetization was computed 
and in a boundary state was found. The system flows from free (with respect to the spin 
variables) to fixed boundary condition that corresponds to boundary spins directed along 
the magnetic field. We show that the boundary entropy monotonically decreases along the 
flow between the fixed points values. Mathematically both models being Gaussian are close 
relatives that is revealed by similar looking answers. However the Ising model besides being a 
more physical example also provides a better illustration of various contact terms appearing 
throughout the proof. These contact terms are reviewed along with the proof of the gradient 
formula in the next section. For the Ising model we check that the gradient formula works 
and demonstrate by explicit computations that it does it according to the general proof. 
To finish the overview of the paper let us mention that appendix A provides a discussion 
of a certain type of distributions supported at a point on the boundary as well as explicit 
computations of the ones that appear in the models considered. Appendix B contains some 
integrals and series used in the computations. 
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2 Review of the proof 

In this section we are going to review the proof of g-theorem that was found in The 
stress-energy tensor contains a boundary piece 

T»v = T b ^ lk (x, r) + 5(ixx)8„ t 5 ut 6(t) . 

Here 6{r) is a boundary operator of canonical dimension 1. The conservation equations read 

&»r™ k {x,T) = o, (5) 

^ u T lk (0,T) + d T 6(r) = 0. (6) 

The trace of the stress-tensor is 

r; = 8 Wk (i,r)+%)e(r). (7) 

Since the system is critical in the bulk Obuik{%, T ) = up to contact terms. The boundary 
trace 9 can be decomposed into a linear combination of the boundary fields a : 

9 = (3 a (\)0 a + h(\)l. 

The generator of dilatations is Tjf so that the renormalization group equation for In z is 

IA—ft— = / / fi 2 dTdx(S(/ix)e(T) + Q bulk (x, t)) = (3 a d a In z + ^h{\) (8) 



dpi 

while for one-point functions of boundary operators we have 

H-fi-(<Pa(n)) = J J li 2 dTdx{[5(iJ,x)6(T) + 6 Mfc (x,r)]0 a (ri)) c = 

P b d b (Mri)) + (l b a -5 b a )(Mri)) (9) 

where 7^ is the anomalous dimensions matrix. 

To prove the gradient formula (jSJ), ® we start with the expression on the left hand side 

a 

g ab (3 b = [ ^dn [ fidr(Mri)0(r)} c (l - cos[27r(r - n )/0\) . (10) 



Since we assume that the UV behavior is governed by some fixed point the singularity in the 
two point function (0 a (ri)^(r)) c cannot be stronger than \t — Ti| -2 . Therefore the integral 
is convergent due to the presence of the 1 — cos[27r(r — Ti)//3] factor. For the purposes of the 
proof we need to split the above expression in two terms: 



P J dr(Mri)0(r)) c (11) 

and 

Aa=-/3 j dr(0 o (T 1 )0(r)) c cos[27r(T-T 1 )//3] . (12) 

At this point the two point function has to be treated as a distribution and the integrals 
are finite due to the presence of contact terms (see e.g. [Ej, section 1.3 for a general 
mathematical discussion) . 
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Integrating by parts 2 in (J12j) we obtain 



A a =^- I fidT(Mri)d T e(r)) c sm 2n(r - Tl )//3 



2tt 

By conservation law (|BJ) 

A a = -2f3 J /x 2 dr(0 o (ri)T^*(O,r)) c i; r (O,r) 

where 

t» T (0, r) = 7- sin [2tt(t - n)//3 

47T L 

is a vector field on the boundary. This boundary vector field can be extended to a conformal 
Killing vector field on the whole half-cylinder. That is 

dpV v + dyV^ = g^daV 01 , v x (0, r) = . 

This vector field generates a subgroup of SL(2, M) -transformations that leave the boundary 
point r = T\ intact. Its explicit form can be most easily found in complex coordinates 
w = 2ti(x + ir) I 'ft. The analytic component reads 

2tt 1 

— (v x + i v T) = _( e — i _ e -»+^ 
ft 4 

and the antianalytic one is given by the complex conjugated expression. The divergence of 
this vector field reads 

d a v a = cos[2tt(Y - n)//3] cosh(2vrx//3) . 

The fact that d a v a (0,Ti) = 1 means that the vector field acts locally around the boundary 
point r = T\ as a dilatation 3 . Integrating by parts in the balk we obtain by (J3J 



A a = 2 J J fx 2 drdx(UnK: lk (x,T)) c d^ . 



There is no boundary term from infinity because of the boundary condition Since t> M is 
a conformal Killing vector we have 



A a = J J /j, drdx(0 a (ri)e bu i fc (x,r)) c (9 a w a . 

Now note that the operator can have only a contact term OPE with any other operator. 
(We are going to discuss these contact terms in more details after we finish reviewing the 
proof.) Since Qbuik has canonical dimension 2 and in a renormalizable QFT <j) a has dimension 
less or equal than 1 the most singular contact terms possible are of the form 8{x)8'{r — T\) 
and 8'(x)S(r — r'). Noting that the difference d a v a — 1 vanishes to the second order at x = 0, 
t = 7~i we see that contact terms make no contribution. Therefore 

A a = // 2 drdx(0 a (T 1 )6 6u ^(x,r)) c . 



2 Since the correlation function is a distribution integration by parts is a valid operation. 

3 Note that there is a one-parameter family of SL(2, M) vector fields that leave a given boundary point intact and act as a 
dilatation in a neighborhood of that point. The particular vector field we chose is the best approximation of local dilatation 
in the sense that d a v a — 1 = 0((r — r') 2 ) that is crucial for the proof. For all other aforementioned vector fields this is just 
O(t-t'). 
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Plugging this back into (JIUj) we obtain 

9abP b = J fJdn J J fi 2 dTdx([S(fix)9(T) + Q bulk (x,r)](p a ( y T 1 )) c 

that by RG equation Q can be written as 

f d d d 

9abP b = J ndrin— {<f>a{n)) = (>— - l)d a \nz = dail^-Q^ - l)\nz = -d a s 

that completes the proof of the gradient formula. 

In the course of proof we encountered contact terms of two kinds. The first kind has to do 
with nonintegrable singularities in the two point functions (0(t)#(t')). Contact terms that 
renormalize these nonintegrable singularities and promote the two-point functions to distri- 
butions come about when in the course of proof of the gradient formula we split expression 
(HOI) into 2 terms (fTT jl . (|T2 |l . 

A textbook example of treatment of such contact terms by standard renormalization 
technique arises when one considers a correlator of composite fields in free field theory 

(:0 2 : (p) :0 2 : (0)) = J d d xt ; »*(:,r: (x) :0 2 : (0)) 

here <p{p) stands for the field momentum space modes. The corresponding momentum space 
Feynman diagram is 




It has a divergence that is essentially of the same kind as a one-loop divergence in an in- 
teracting theory. The usual momentum space counterterms when translated into position 
space correspond to a linear combination of derivatives of a delta function with divergent 
coefficients. (See e.g. section 6.2.2 for a sample computation.) Alternatively instead 
of working in momentum space one can employ a point splitting regularization + minimal 
subtraction to compute these contact terms. In the free massless boson with a boundary 
mass model the two point functions are integrable and there is no need for contact terms 
of this kind. In the Ising model with a boundary magnetic field the two point function has 
a nonintegrable singularity ~ |r — Using the point splitting plus minimal subtrac- 

tion scheme the distributional two-point function (0(r)^(r')) is computed in section 4 (see 
formulas ijSBI) . (153)1 ). 

Note that although such contact terms are scheme dependent the sum of terms (fTTj). (fT2|) 
is independent of the choice of contact terms. 

A contact term of a different kind is a contact term in the OPE between Obuik{%, T ) and a 
boundary operator 0(r'). In a bulk theory analogous contact terms stem from the fact that 
®buik{z, z) is a generator of dilatations. For instance if Q(z, z) is a CFT primary of weights 
(A, A) we have a contact term 

Q hulk (z, z)®{z\ z') = -(A + A)5^ (z - z')$(z', z') . (13) 

One can derive this contact term starting with standard OPE of $(2, z) with T(z) and T(z) 
and using the Ward identity 

2[dT zz + dT- zz }(z,zMz',z') = -5^{z-z')d^{z',z!) 
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and the complex conjugated one. Namely if we differentiate the OPE of T(z) and $ with 
respect to 8 and make use of the identity 

8—^— = it5 {2 \z-z') (14) 



we obtain 



dT zz (z)<S>(z',z') = ^-d5 {2 \z - z')$(z\z') - \b {2 \z - z')d$(z',z') 



' z') 



Plugging this into the above Ward identity we deduce that 

2T zz (z,z)®(z',z') = -A5<- 2 \z- z')$(z 

A complex conjugated identity is obtained in an analogous way. Since 6 = 2{T ZZ + T zz ) we 
obtain (JT3J). 

In the presence of a boundary the generator of dilatations becomes Bbulk( x , T ) +S(x)8(t). 
The Ward identity corresponding to translations preserving the boundary reads 

[*-(d - 8)& bulk + 2t(8T zz - dT- Z z)W) + 6(x)[i(T zz - T- zz ) + d T 9](f)(d') = 
5{x)5(t - r')d Tl (j){r') . (15) 

We see that there is a freedom in how the contact term on the right hand can be split between 
the bulk and boundary 4 terms on the left hand side. In other words the contact terms between 
(f) and bulk and boundary counterparts of the stress tensor are scheme dependent. Given 
a particular renormalization scheme we can again take OPE's of T(z) and T(z) with 0(r') 
and differentiate them with respect to z and z respectively using identities similar to ()14j) . 
Since for any z ^ it' the OPE's with T(z) (T(z)) are holomorphic (antiholomorphic) in z (z) 
a derivative with respect to z (z) taken in a distributional sense should yield distributions 
supported at a boundary point r'. A more detailed discussion of how to obtain formulas 
similar to (J 14)) in the presence of boundary is given in appendix A. Once we have a contact 
term of with dT and dT we can use them in (|15)l to obtain a contact term of Qbuik with 
<p. We follow this sketch of a computation for the two models we consider in the subsequent 
sections. Those computations are done in a particular renormalization scheme and we would 
like to emphasize again the scheme dependence of contact terms. Thus one could make use 
of a scheme in which there are no contact terms in the Qbuik OPE with <fi at all with all 
proper contact terms appearing in the OPE between 9 and <j). 

Note that contact terms of both kinds are not independent. They appear together in the 
RG equation for a one-point function of 

V^rW)) = I fidr{<P{T')e{T)) c + [ [ /i 2 rfa;dr(0(r')e^(x,r)) c . (16) 



<9/i 

When the subtraction scheme is chosen one can compute the left hand side. If the scheme 
for computing the contact terms of both kinds in the correlators is chosen consistently then 
upon integration (jlfij) should hold. Another appearance of both contact terms happens in 
the identity 

/^t(0(t')(9(t)) c cos — (t-t') (17) 
that was derived in the course of proving the gradient formula. 

4 By a boundary term here we mean the term proportional to S(x). 
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3 Free boson with a boundary mass term 
3.1 The gradient formula 

As a warm-up example we are going to consider in this section a free boson with a boundary 
mass term. This model was first considered in 0. The action functional on an infinite 
half-cylinder reads 

P oo P 

S = — [drfdx <9„0<9> + — I ixdrcf) 2 . (18) 
8vr J J 8tt J 

oo o 

Here \i is a renormalization scale inserted so that the boundary coupling constant u is 

dimensionless. Varying this action functional we obtain a boundary condition known in 

mathematical literature as Robin boundary condition: 

-|^(t,0)+«/^(t,0) = 0. (19) 

The field theory space for this model is labelled by a single coordinate u conjugated to a 
boundary operator 



1 ..2. 



Noting that 



87T 

T 

J- 7.9 



47T 

it is easy to derive using the boundary condition ()19j) 

-l 



(T zz -lk)| x=0 = -i^—d T ■ (20) 

Therefore by © the connected part of the operator 9 is 

u 

C = ~^-4> 2 = Pu4> u 

where (3 U = u is the beta function of the coupling constant u. 
The disk partition function was found in [S] to be 

z = i v /^e^r(a)e- ffln ^/ 2 ^ (21) 

where a = ufi/3/2ii and T, 7 are Euler's gamma function and constant respectively. The last 
factor in the above expression is absent in Witten's computations due to the specific choice 
of the renormalization scale /i. Also we included a normalization factor l/2y / 7r that can be 
computed by a method similar to the one used in '13,.. From the RG equation 



d In z dlnz 
^ = U ~du~ + W 



we find the whole operator 9: 



e = 9 c + hl = -^: ( p 2 :-^l. (22) 

877 ZTC 



5 The unusually looking factor of fi 2 in this equation is due to our conventions in which all operators are dimensionless. 
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From (|21j) the boundary entropy is readily found to be 

s = - ln(2v^F) + - ln(o-) - - + In T{a) - <np(a) + a . (23) 

The boundary entropy (}2*3*j) was first computed in jU] and its monotone decrease with a was 
checked. In the limit a — > oo the boundary entropy s — > sjr = — ln\/2 that is the value 
corresponding to Dirichlet boundary condition. 
We want to check the gradient formula 

a _ 
Quuyu 

ou 

Here 

9uu = J [idn J fidr 2 (l - cos[27r(r 1 - r 2 )//3])(i-</» 2 (r 1 )i-</» 2 (r 2 )) c (24) 
o o 

where the correlator is normalized and connected. The two point function was computed in 
On a half-cylinder it reads 

((j){zi)<t>{z 2 )) c = - In \ e - 2nZl/ ? - e - 2nZ2/ P\ 2 - In 1 1 - e -M*+*Vf>\ + 

O °° p -27rn(2i+2 2 )//3 , -27rn(2 1 +2 2 )//3 

• ^ 

a ^-^ n(n + a) 

n=l v ' 

The short- distance divergence in the correlator is logarithmic and thus integrable. We can 
then split the computation of ()24j) into computing separately the term with identity and the 
term with cosine. From (|2*3j) we obtain 

^| rfri | rfr2COs[ 2^ (ri _ r2 )/^ ( ^ (ri) ^ 2(r2))c= h*£j I ; (26) 



^/m£*"W^M>.= g^-^ (27) 



Combining (J1S1) and (J2*7)l we obtain 



2-k J L 2a 

that using (|2*3j) can be readily checked to coincide with —d u s. 



3.2 Contact term in the OPE of <p u with Qbuik 

It came out in the above computation that the term (J26|) coming from integration with cosine 
in the metric appears on the right hand side of the gradient formula from the In \i term in the 
exponent of (|2ip. This term comes from normal ordering subtraction defining the operator 
■.cf) 2 :. Another way to look at this logarithmic divergence is that it arises from mixing of 
operator : cf) 2 : with the identity operator under a scale transformation. Indeed according to 
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the proof of g-theorem reviewed in section 2 the term (|26|) should reproduce the bulk part 
of the dilatation operator acting on <f> u . As it was already mentioned in section 2 this term 
can be fixed by the RG equation ([160. Assuming that <p u has no contact terms in its OPE 
with 9(t) 6 we can compute the first term on the right hand side of (fT""j). it is given by (|27jl . 
The left hand side of (fTf"") can be computed using (|2*Tj) as 

It can be checked then that the integrated OPE with Qbuik does match with the cosine term 
in the gradient formula that is (f""""|) multiplied by a factor of u. It follows that the OPE itself 
is fixed assuming that it takes up the form 

Q bu ik{z, z)4> u {t') = Const 5 {2 \z - ir')l . 

It is instructive to derive this contact term directly in a certain regularization scheme 
using local conservation equations. Consider the bulk-to-boundary OPE of T zz (z) with the 
boundary operator (j) u . We are interested in the most singular part of the short distance 
expansion. It can be computed using the bulk-to-boundary propagator with the Neumann 
boundary condition. We have 

T sg = — : d<pd<i) : 
An 

and the OPE 

T zz (z)<f> u = --^ :«90<90(z): :0 2 (r'): = 

™ V 1 + ., * fg W * , + - (*» 



167r 2 V 1 - e- 2 < z ~ iT ')l^ J 27r 2 V (3 J z 1 - e- 27r ^" ir ')/' 3 

where dots stand for less singular terms. We further differentiate this term with respect to 
z using the identity (jA.9|) 



1 _ e -2nz/P y 2 7f J 2 

that is derived in appendix A. We obtain thus 

BT zz {z)(~) : ( f?:{T') = ^-d z 5 {2 \z-iT')l + ... (29) 

and the analogous complex conjugated identity. 

Substituting ([2*9*)) and the complex conjugated expression into the bulk part of (fi3j) we 
obtain the OPE 7 

e b uik(z,z)Mr') = Qbuik(z,z) (--M :</» 2 (r'):= -^ 2 \z - ir')l . (30) 

Using this contact term and formulas (|2*7)l . (|2*B|l . (|2*T|) we find that the identities (fTHj) and 
fjTTf) hold. 

6 This should be considered as part of the particular renormalization scheme we choose. 

7 Note that the less singular terms in 1291 proportional to the identity operator drop out when we take the difference 
9T zz <f> u — dTzz4> u while the terms proportional to 8 T (j> u coming from single contraction terms in i'liii combine with analogous 
term in the boundary-boundary OPE in 1151 to yield the right hand side. 
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4 Ising model at criticality with a boundary magnetic field 



4.1 The model 

A more physical and also more illustrative example of a boundary flow is a 2D Ising model 
at criticality perturbed by a boundary magnetic field. The system flows from a free spin 
boundary condition in the UV to the fixed spin boundary condition in the IR. On an infinite 
half-cylinder the action functional of the model reads [11. 

S = — dxdT(tfjdtfj + tpdip) + / dT(^-ifjtfj + -aa + ihfi 1/2 a(ujtfj + uii!)). (31) 
2n J J J An 2 

Here lu = e^ 4 , u> = e~ J7r//4 , the complex coordinate is z = x + ir. The boundary fermion 
a(r) accounts for the double degeneracy of the ground state. Its free propagator reads 

(a(T)a(r')) free = ^sign(r - r') . 

The boundary coupling constant h is dimensionless. The equations of motion and the bound- 
ary conditions read 

dip = , dip = (32) 
a = — ihfi 1 ^ 2 (ujip + uji/j) , (33) 
uip - ui> = -Anhfx 1/2 a . (34) 

Functionally integrating out the boundary fermion a we obtain a functional integral over 
the ip, ip fields of the form 

D[a]D[i),ij\e~ SbM ~ Sd = 
J D[V>,$]e- Si »»*exp J j drdr'^(T)sgn(T - r , )*(r / )^ (35) 

where 

$ = Dip + OJIp . 

The boundary term in (|33j) imposes a local boundary condition on the remaining fermionic 
fields 

(d T - i\)ujip = (d T + iX)^ (36) 

where 

A = 4nh 2 fi . 

Note that on a finite length cylinder two sectors: Neveu-Schwarz and Ramond are present 
in the bulk theory. In the infinite cylinder length limit L — > oo the Ramond sector states are 
suppressed by a factor exp(— L/16/3) due to the Ramond vacuum weight of 1/16 and thus 
these states do not survive in the thermodynamic limit we are considering. We will thus 
consider below only the Neveu-Schwarz sector. 

It is easy to find the two-point free fermion Green's functions satisfying the boundary 
condition 



<V(*lW(*2)> 



sinh \j${zi - z 2 )j 
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WfaMh)) = j^- v + g(z 1 + z 2 ) (37) 

smh \^{z\ + z 2 )j 

where 

4.2 The operator # 

The connected part of the boundary operator 9{r) can be found from conservation equation 
(0). We have 

T„(z) = -^T(z) , T(z) = 

where the unusual factor of fi~ 2 is due to our scaling conventions. Thus on the boundary we 
have 

T XT (0, r) = i{T zz - T 2 z)\x=o = ^(ipdrip + ipd T iP) . (39) 

The boundary condition (|3l)Jl then implies 

d T ip = - i(d T + i\)i> , (40) 
d T tp = i\$ + i(d T - i\)i> . (41) 

Using these formulas we obtain 

T XT = ~d T {H>) (42) 

and therefore 

C =^:#:. (43) 

47T 

We will derive shortly the term in 9 proportional to the identity operator. 

Another way to find 9 C is by applying Noether's theorem to the Lagrangian (|31|). Taking 
the classical mass dimension of if) to be 1/2 we have 

9 C = — ia{ujip + uip) . 

On the boundary condition (|36|) this expression coincides with (|43|) . The last expression 
however manifestly has the form f3h4>h{ T ) with flh = h/2 and 

<p h = -iafi' 1 (utfj + uii/j) = : -0$ : . 

4.3 The partition function 

The disk partition function for the model at hand was found in ()A.7|) using the boundary 
state formalism. To fix the nonuniversal term present in the partition function in a particular 
subtraction scheme (to be below) we derive the partition function using the correlators (J37|) . 
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It follows from ()35|) that up to a normalization constant the disk partition function can be 
extracted from the equation 



d\nz(h 2 jj,) 1 



dh 2 n 



~j J c/rrfr , (^(r)^(r / ))sgn(r-r / ). 



(44) 



Using ()3Tj) the two-point function in the last expression can be computed to be 

(*(t)*(t')> = i[<V(r)V(0) - $(tMt'))} + (^(r)^(r')) + (^(r)^(r , )> = 
4vr//5 



sm 



E 

n=0 



4A 

n + 1/2 + /3A/2vr 



sm 



^(n + l/2)(r-r') 



(45) 



The contribution of the series to (|44j) can be expressed via Euler's psi function while the first 
term contributes an integral 



1^- 



(3 P 

dr / dr'- 



o 



o 



sm 



(46) 



that is logarithmically divergent. As shown in appendix A (formula ([A. 15)) ) the point split- 
ting regularization plus minimal subtraction yield a renormalized value 



/?(21m>/3)-V>(l/2)). 



(47) 



Collecting all terms we obtain 
din z{h 2 ) 



dh 2 



/?(21nGu/?)-V>(l/2)) + X; 



Ah 2 j3 2 n 



n=0 



n + l/2 + 2(3h 2 /j,)(n + l/2) 



(32 ln(/3/i) + Pi>(l/2) - ip{2(3h 2 (i + 1/2) . 



This yields 



C 



D a(ln(/3^)+|V(l/2)) 



(48) 

(49) 



r(a + 1/2) 
where 

a = 2(3jj,h 2 , 

C is a normalization constant. Up to this constant and nonuniversal terms of the form e aC> 
that depend on the subtraction scheme the above value of z coincides with the one computed 
in The constant C can be fixed by computing a partition function on a finite cylinder 
that can be canonically normalized and by studying its factorization in the infinite cylinder 
limit. This yields C = y/rc [13 . 

The logarithm of partition function (|49|) satisfies the RG equation 



/'- 



d In z hdlnz 



dfx 



2 dh 



+ 2/3/i/i 2 . 



Thus the complete 6 reads 



6 = —ii- 1 :^: +2(3rih 2 l 

47T 



(50) 



(51) 
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From (J49|) one readily computes the boundary entropy 

s(a) = lnv^F-lnr(a + l/2) + a^(a + l/2)-a. (52) 

One finds then that suv — s (0) = 0, sir = s(oo) = — \ In 2 in accordance with the boundary 
entropy values for free and fixed boundary conditions in the Ising model jTOj. 



4.4 (0/ l (r)6 , (r')) c correlator. A contact term. 

We are interested in this section in the correlator (4>h{T)0(T')) c and particularly in the contact 
term it contains. Up to a constant factor this correlator equals 

(:^: (r) :#: (r')) c = G(t - r') . (53) 

Let us write 

7T//3 



(V(r)V(r')) 



sin 



+ ^(r-r')] 



(54) 



where the function g is defined in (jHHjl . Then by Wick's theorem 



G(T-T / ) = 0/[<(r-r , )]-ff[<(T / -r)])- 



7T//3 



sin 



g[i(r - r')\g[i{r' - r)\ . 



It is easy to check that the second term is an integrable function while the first term contains 
a singularity of the order |r — r'| _1 . This is unlike in the boundary mass model where the 
correlator (8(t)4> u (t')) c was integrable and there was no need for a contact term that would 
ensure the integrability. After some work we can single out the singular piece in a simple 
form and rewrite G{r — t') as 

1 



G(t - t') = -2A 



7T 



sin(^(r-r'))| 



+ fi(r ~ r') + / 2 (r - /) 



where 



fi(r-r') 



2A v ^ 



sin 



E( 2n +1)(T-T') 



sin(f (r - r')) ^ (n + 1/2 + a)(n + 1/2) 



/ 2 (r-r')=4A 2 £ 



cos Rf- (t — r')(m — n) 



(55) 



(56) 



(57) 



n (ri+l/2 + a)(m + l/2 + a) ' 

n,m=0 v ' ' v ' ' 

Both functions f\ and fi are integrable. Applying a point splitting plus minimal subtraction 
to the singular piece (see formula ()A.15|) in appendix A) we obtain a distribution that we 
denote 

R 1 

f|sin(|(r-rO)| 

times a factor — 2A7r. This distribution acts on a test function 4>(t) periodic with period (3 

as 



(R 



1 



^|sin(|(r-r'))r r/ 

0(r) + 0(/3 - r) - 20(0) 



sin(|(r-r')) 



+ 0(O)(21n(^)-^(-)). 



(58) 
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Thus we promoted the correlation function G(t — t') defined for finite separations in (J55I) . 
f)56|) . (J57)l to a distribution 

G(t - r') = -{2\k)R , 1 — + A(r - r') + / 2 (r - r') . (59) 

§|sin(|(r-r'))| 

Now we can compute the integrals involving the distribution (JSTJj) . A straightforward 
computation yields 

r-/3 



jf J drdr'G(r - r') = 16vr 2 fy,( a + 1/2) - ^(1) - | ln(/5/i) + a^'(« + 1/2)) 



/3 /■/? 




drdr'G{T — r') cos 



Jo 



2tt 
7" 



r — r 



1 , ,1. 



167r 2 a l-i/,{a + 1/2) - -^(-) - - ln(/?/i) + 1 

Note that the renormalized value of the integral ()60|) coincides with the value of 

2 2 d 2 Inz 



(60) 



(61) 



-2niYh 



dhdh 



This of course was to be expected because we used the same point splitting + minimal 
subtraction scheme when evaluating both quantities. 



4.5 The gradient formula 

We would like to check now the gradient formula 

ds 



dTdr'((j) h (T)6(T')) c (l - cos 



o o 



2vr 
J 



[T — T 



(62) 



The RG equation implies that 



ds 2 ds 2/r d In z 

dh h^dfi h dfidfi 
In the last subsection we computed 

2/x- 1 1 



: (:^:(r):^:(r')> c 



2/i- 2 1 



T — T 



h 10~ 2 "" h 16~ 2 

where G is a distribution given in (|59|). Checking the gradient formula (jo^j) then boils down 
to checking 

fi 



, d 2 In z 



<9/t 2 



16tt 2 



drdr'Gir — r')(l — cos 



o o 



2vr 



r — r 



(63) 



The left hand side is immediately computed using ()49j) with the result 

n 2 



dfx 2 

that indeed equals to the difference of integrals (|5U|) and (|BT|l times a factor of —(16 



2\-l 
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4.6 Contact term in the OPE of 6 h with 



bulk 



As in the bosonic model we can derive the contact term in the OPE of Obuik with <ph by 
using Ward identity (Ti3j) and the OPE's of <ph with T(z) and T(z). The term proportional 
to the identity operator in the OPE of T(z) with iipip: (r') reads 

0) (?) = l[^{z)^{T'))d^{z)^{r')) - (^(z^^d^z)^))] + ... 

The most singular terms of the order (z — ir)~ 3 in the above expression cancel out and the 
leading singularity comes from the logarithmic subdivergence in (ipip). Up to less singular 
terms we have 

:m: (*) *fr (f) - 1 (f ) ft^&ff +.~ (64) 
As z — > ir the function g(z — it') has asymptotics of the form 

g ( z - ir') = 2i\ ln(l - e -M*-^)/P) + 2 i\[ 1 + ^{a + 1/2)] + 0{z - ir') 

where 7 is the Euler's constant. Differentiating the above OPE with respect to z in the 
distributional sense we find using formulas fO|l . (|A~T6jl . (IB~4|l that 



:ipdip: (z) \ ipip: (r') 



l(27riA) [^(« + 1/2) - ^(1/2) - 1 ln(///3) + l]d z 5^(z - zr') + . . . (65) 

where the dots stand for less singular terms proportional to 8{z — iT')d<ph- We also have an 
analogous contact term with dT. Using these contact terms and the Ward identity (fT3j) we 
obtain a contact term with Qbuik of the form 

<S>buik(x,r) -.ip'tp: (t') = 

l(-8rf)[^(« + 1/2) - ^(1/2) - \ ln(///3) + l]5(x)5(r - r') . (66) 

Using this expression along with (jBHj) . (jBTj) we find that equations (fTBj) . (|T7j) indeed hold. 
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A Some distributions supported at a point on the boundary 

In this appendix we discuss how one can define distributions formally expressed as df(z) 
where f(z) is a function holomorphic in the x > region of the cylinder and contains a 
singular point on the boundary that without loss of generality can be chosen to be x = 
0,r = 0. We further assume that f(z) diverges at that point slower than z~ 2 . In particular 
we will be interested in defining 

8 T^h>? (A - 1) 
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and 

ln(l-e-^) 
^ l_ e -2W/3 ' [A - Z) 

It is intuitively clear that any distribution expressed as df where / satisfies the above 
conditions should be supported at the singularity point z = 0. By dimensional reasons the 
answer should be proportional to 8^ 2 '(z). To obtain a rigorous definition one can proceed 
as follows. Let <p{z, z) be a test function which is smooth and fast decreasing at infinity. 
Consider a formal identity 

oo P oo/3 P P 

Jdxjdr (pdf = - Jdxjdr (pdf + ^ Jdrcpf\ x=00 - ^ Jdripf\ x=0 (A. 3) 

00 00 

that follows from a formal application of Stokes formula. The first term on the right hand 
side is well defined as the singularity of / is integrable in two- dimensions. The second term in 
the above equation vanishes because the test function is fast decreasing at infinity. In the last 
term one has to regularize a singular function of one variable f(x = 0, r) and promote it to a 
one-dimensional distribution. Thus we see that the problem of defining a distribution df on 
a space with boundary (a half-cylinder in our case) boils down to defining a one-dimensional 
distribution to be denoted Rf(r) that regularizes the restriction of / to the boundary. This 
can be done in a standard fashion, see e.g. [Hj. As we stick to point splitting regularization 
throughout the paper we will define a distribution Rf(r) by the formula 

P P-e 

JdrcpRf := lim[ J dr V f\ x=0 - C(e)(p(0)) (A.4) 

0+e 

where (7(e) is a counterterm ensuring the convergence of the whole expression. In the case 
when no counterterm is needed and it can thus be set to zero the right hand side gives the 
principal value of the corresponding integral. When (7(e) is present it is defined up to an 
essential ambiguity (7(e) — > (7(e) + Const. 

As the first term in (|A.3|) is well defined we can rewrite it using Stokes formula as 



00 P P-e 

/dx drcpdf = lim / / dxdnpdf = — lim[ / dr(pf\ x=0 + / dzcpf] . 

J e -*°J J\z\>e ' 2 e^O J J\z\=e,x>0 

0+e 

Substituting the last expression and expression (|A.4|) into (|A.3j) we obtain 

00 p 

fdxfdripdf := Jlim[C(e)</?(0) + / dzcpf]. (A.5) 

J J 2 J\z\=e,x>0 

11- 

This expression makes it manifest that 

Bf = C f 5 {2 \z) (A.6) 

and that up to the already noted ambiguity in (7(e) the coefficient Cf depends only on the 
behavior of / in the vicinity of z = 0. For practical purposes however it is convenient to use 
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a different representation of the constant Cf. By applying Cauchy formula to the function 
/ in region \z\ > e we have 

p-e P 

I dz f = - drf\ x=0 + drf\ x=O0 . 

J\z\=e,x>0 J J 

e 

Using this formula we deduce from (|A.5|) and f)A.4|) that 

C/ = -i JdrRf + ^ jdrf\ x=OD . (A.7) 
o o 

Let us apply now this formula to functions (jA.l|) . (jA.2j) . The restriction of function / in 
(jA.lj) to the boundary reads 



1 1 % ( 



TXT 



Ur) = .^-^ = cot — . (A.8) 

1 - e-W/i 2 2 V P J 

It is clear from this expression that the principal value integral of /i(t) exists. We thus set 
= P.V./i. Taking into account a nontrivial contribution at infinity one readily obtains 
from (TP)) 

8 = -fi®(z) . (A.9) 

This formula is to be compared with standard formula derived in the absence of boundary 
when z = z' is a point in the bulk: 

d- z —^— = Tx5 {2 \z-z'). 
z — z' 



We see that in our regularization scheme (jA.9|) corresponds to having a half of delta function 
in the case when the singularity point is on the boundary. 

To derive a distribution corresponding to (jA.2|) we first note the formula 



ln(l - e- 2 ™ 115 ) = ln(2| sin(vrr//3)|) + % -{tx - 2nr/p) (A. 10) 

where it is assumed that < r < (3. We further notice that the function / entering (|A.2|) 
vanishes at infinity. Thus to find the corresponding C/ it suffices to define a regularized 
integral 

p 

ln(l - P" 27r W/3^ 



Cf = — dr R . 

' 2 J 1 _ e ~2mr/p 



Multiplying (|A.8jl by ()A.10jl and dropping the terms vanishing under the principal value 
integration we arrive at the expression 



p p 

dr ln(2| sin(7rr//3)|) - - [drRcot ( ^ ] (2ttt//3 - tt) 



27 v/?7 
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Keeping the limits of integration regularized we can rearrange the terms separating a finite 
and a singular piece in this integral as 



p/2 
dr 

0/2 

71 



;^ms 7TT" . . „. /7T\ COs(71t/3) — 1 

-In 2 sm 7TT//3 + — cot ttt//3 + - , l , 1 P 
1 2p \4/ sm(7rr/p) 



dr 



4 J sin(7rr//3) ' 



(A.11) 



The integral on the first line converges in the limit e — > and can be evaluated using formula 
f)B.2j) . Namely the value of the convergent term in (jA.llJ) equals (/3/2) In 2. We can extract 
a finite value out of the divergent term 



p/2 



dr 



sin(7rr/ (3) 



(A.12) 



by adding a minimal counterterm: 



C(e) = ^(ln(e)+lm>/3)) 



7T 



where /x is the renormalization scale. The remaining finite piece can be computed using the 
following trick. Consider an integral 



do 



e 

sm" | - 

2 



V > — . 
2 



(A.13) 



The minimally subtracted value of this integral at v = —1/2 can be obtained by taking 
v = — 1/2 + e, e > 0, expanding in powers of e and subtracting the pole. This yields the 
value 

J 2 (-l/2) = -V 1 ' 



where ^ is the logarithmic derivative of the Euler's gamma function. This fixes the finite 
part remaining after the minimal subtraction in the divergent term (|A.12|) . We have 



13/2 

JdrR 

o 



P/2 



lim 



sin(7rr//5) e-*oLj sin(7rr//3) 



C(e) 



A)[21 n (^)-^(l/2)]. (A.14) 



In a similar fashion a minimally subtracted value of the following integral that is used in the 
main body of the paper is obtained 



Iren — 



dr I dr'R- 



1 



o o 



sm 



r — r 



[21n(/i/3)-^(l/2)]. 



(A.15) 
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Collecting all pieces together in (jA.llj) we finally obtain 

9 H i~^ZT = (I) -h 2 -\ in( ^ } + ^ (i/2) ] 6(2) {z) • (A - i6) 

Note also that applying the general formula (jA.7|) we find 

<91n(l - e - 27Tz/l3 ) = 
that agrees with the naive dimensional intuition. 

B Some useful integrals and series 

sin(0/2) V ' 



dx ln(sin(x)) = -tt In 2 . (B.2) 



x i 

ip{x) = - 7 + (x-l) Y— -. (B.3) 

v ' v ' ^ nix + n-l) 

n=l ' 

^(1/2) = -7 -2 In 2. (B.4) 
V Sm((2n + ; 1)T) = tt/2 for < r < tt , -tt/2 for tt < r < 2vr . (B.5) 

n=0 ' 
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